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Abstract. We study the existence and uniqueness of nontrivial stationary solutions 
to a nonlocal aggregation equation with quadratic diffusion arising in many contexts in 
population dynamics. The equation is the Wasserstein gradient flow generated by the 
energy E, which is the sum of a quadratic free energy and the interaction energy. The 
interaction kernel is taken radial and attractive, nonnegative and integrable, with further 
technical smoothness assumptions. The existence vs. nonexistence of such solutions is 
ruled by a threshold phenomenon, namely nontrivial steady states exist if and only if 
the diffusivity constant is strictly smaller than the total mass of the interaction kernel. 
In the one dimensional case we prove that steady states are unique up to translations 
and mass constraint. The strategy is based on a strong version of the Krein-Rutman 
theorem. The steady states are symmetric with respect to their center of mass xq, 
compactly supported on sets of the form [xq — L, xq + L], on their support, strictly 
decreasing on {xo,xo + L). Moreover, they are global minimizers of the energy functional 
E. The results are complemented by numerical simulations. 



1. Introduction 

Phenomena with long-range aggregation and short-range repulsion arise in many in- 
stances in population biology such as chemotaxis of cells, swarming or flocking of animals. 
A variety of mathematical models has been proposed for such situations, at the particle 
as well as at the continuum (mean field) level. In particular, if nonlocal repulsion acts at 
a smaller scale with respect to nonlocal attractive forces in the large particle limit, then 
a nonlocal repulsion term can be replaced by a local term with nonlinear diffusion, we 
refer toIllllIIlEIllMllSnilSSlEIlliniliniEelfor several examples. A prototype model, 
which we shall also investigate further in this paper, is given by 

dtp = div{pV{ep-G*p)) (1.1) 

where the convolution is carried out with an aggregation kernel G such that G{x) = g{\x\) 
with g'{r) > as r > 0. This model arises in a natural way as the limit of a stochastic 
interacting particle model with pair interactions (cf. \43\ H3] respectively |26j for general 
background). Models with the same structure have been recently used to model opinion 



formation, cf. \48\ I50j . Here we shall restrict to the case in which (1.1) is posed on the 
whole space W^. 



In case G is A-convex, then the equation (1.1) can be formulated as a gradient flow in 



the Wasserstein metric (cf. [H [SH [55j ) of the energy (or entropy) functional 

^ / p^{x)dx - J 
^ ./rod ^ 



E[p]:=^ [ p\x)dx-\l / G{x-y)p{y)p{x)dydx, (1.2) 
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see [m [3l [T71 |T8]. Models of the above form have been investigated with respect to 
several aspects, e.g. existence and uniqueness in the context of entropy solutions [HI [8], 
well-posedness in the context of Wasserstein gradient flows as a special case of the theory 
developed in yj, and - in particular in connection with the classical Patlak-Keller-Segel 
model for chemotaxis |471 132] - with respect to blow-up vs. large-time existence, cf. e. 
g. \30\ [271 m [21] for models with linear diffusion and |33[ [T5] for models with nonlinear 
diffusion. 

An interesting and important question is the characterization of large-time behavior 



of solutions to equations of the form (1.1), which is related with the possible existence 
of nontrivial steady states, even when the quadratic diffusion is replaced by a more 
general nonlinear diffusion. This issue is solved in detail for purely diffusive equations, 
in which solutions decay to zero with a prescribed rate for large times and behave like 
the (compactly supported) Barenblatt profiles, cf. the classical works of Vazquez on the 
self-similar behavior of the porous medium equation, which are nicely collected in the 
book [53], as well as the papers [l5l[T9]. In the purely nonlocal case, namely when e = 0, 
this issue has been studied extensively in many papers [371 [13 133 [3 [H [13 [3 13 [23 123 
[23], combined with the study of the regularity of solutions compared to the attractive 
singularity of the interaction kernel. In particular, solutions are known to concentrate to 
a Dirac delta centered at the initial center of mass (invariant) either in a finite or in an 
infinite time, depending on the properties of the kernel G at x = 0. When the kernel G 



is supported on the whole space, the Dirac delta is the unique steady state of (1.1) with 
unit mass and zero center of mass. 

The asymptotic behavior in the general case with both nonlinear diffusion and nonlo- 
cal interaction has been only partially addressed. A first attempt in this direction was 



performed in [13j, in which the existence of steady states of (1.1) for sufficiently small 
e and the non-existence for large e in the one-dimensional case was proven by means of 
the pseudo-inverse representation of the Wasserstein distance. More refined results in a 
similar model derived in [2g| BB] with cut-off density have been found in [14:\ [M] . Parallel 
to their work, the authors of the present paper recovered the results in [2], in which a 
quasi sharp result of existence of minimizers for the energy E[p] in a multi dimensional 
framework has been proven. 

A key open (to our knowledge) problem in this context is the uniqueness of steady 
states under mass and center of mass constraint, its main difficulty being the fact that 
the functional E is neither convex in the classical sense nor in the displacement convex 
sense [ 41] (except when G is concave on M*^, see [17j ) when e < f G with G > 0. 

In this paper we further investigate the detailed structure of steady states in one 
space dimension. We remark that our work does not go into the direction of providing 
sharp regularity conditions on the kernel G. Roughly speaking, G is smooth, radial 
(with decreasing profile), nonnegative, integrable and supported on the whole space. The 
precise assumptions on G are stated at the beginning of Section [2] We found out that 
for a kernel G decreasing on x > 0, the L^-norm of G compared to e marks a threshold: 

• U e > f G, then there exists no steady state. 

• li e < J G, there exists a stationary state. 

The same results are recovered in [2], except for the critical case e = J G. 

The main result of our paper deals with the case d = 1. Here we can give a detailed 
characterization of the stationary states in the case e < J G. li G' only vanishes at zero. 
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then there exists a unique stationary state (up to translation), which is a minimizer of 
the energy E at fixed mass. The stationary states have compact support, which increases 
with e. Moreover, such a steady state is symmetric with a single maximum at the center 
of mass, and monotone on both sides of the center of mass. The precise statement can be 



found in Theorem 4.13 If e is small enough, the stationary states are concave on their 



support. The main tool in the proof of the main result is the statement of the stationary 



equation as an eigenvalue problem, see formula (4.10), together with the use of a strong 
version of Krein-Rutman theorem (cf. Theorem 4.10), which allows to characterize the 
steady states as eigenfunctions corresponding to a simple eigenvalue. 

The uniqueness of the stationary state is somewhat surprising, since the energy func- 
tional E is not convex for e < J G and thus one might expect other stationary points of 
E. On the other hand also in the case e = one can see that there is a unique (measure) 
stationary state concentrated at the center of mass if G has global support. Adding the 
squared norm for positive e makes the functional closer to convex, and thus probably 
does not lead to additional stationary points. 

The paper is organized as follows. In Section [2] we recall the statement of the problem 
and provide some preliminary regularity results. In Section [S] we complement our results 
with those proven in [2j and provide sharp conditions on e and / G for the existence of 
non trivial steady states. In Section [4] we prove our main results about the uniqueness 
of steady states in one space dimension. Finally, in Section [5] we complement our results 
with some numerical simulations. 

2. Preliminaries 
We consider the evolution equation 

dtp = div{pV {ep - G* p)) (2.1) 

and its stationary version 

= dW{pV{ep-G*p)) (2.2) 

posed on the whole space M*^. Due to the applied setting of the model, we shall consider 
here only nonnegative solutions p > 0. 

Assumptions on the kernel G. We shall assume throughout the paper that the 
interaction kernel G satisfies 

(1) G > 0, and supp(G) = M"^, 

(2) G E W^'^iR'^) n L~(M'^) n C2(M'^), 

(3) G{x) = g{\x\) for ah x £ R'^, 

(4) g'{r) < for ah r > 0, 

(5) g"iO) < 0, 

(6) limr^+oo gir) = 0. 

We emphasize here that providing sharp conditions on the regularity of G is not a 



purpose of the present paper. Let us recall that the equation (2.1) preserves the total 
mass 

M 



and center of mass 



jpix)dx 
CM[p{t)] := J xp{x,t)dx. 
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Moreover, it is easily seen that, for a given stationary state p solving (2.2), Mp and 
p{' + xq) are still stationary states for all M > and xq G M. 

We shall therefore assume M = 1 for simplicity. For future use we introduce the space 



V 



p(x)dx = 1 



Moreover, from now on we shall assume for simplicity 



G{x)dx = 1. 



This is not restrictive since the kernel G can always be normalized by modifying the 
diffusion constant and the time scale as follows 



drp = dw\^pV{e'p-G* p)j , r = G = G/\\G\\Li^^d), e' = e/\\G\\mM.i)- 

Let us recall the following results on the existence and uniqueness of gradient flow 
solutions to (2.1), which follows from the theory developed in [1]. In this sense, let us 
introduce one of the most important tools related with the study of the evolution equation 

e 

p'^{x)dx 



(2.1) and in particular with the stationary version (2.2), namely the energy functional 

e 



E[p] :-- 



G{x - y)p{y)p{x)dydx. 



Theorem 2.1 ([Ij). Let po € nV such that p >0 and E[pq] < +oo. Let G satisfy the 
above assumptions. Then there exists a unique weak solution p to (2.1) with 

• E[p{t)] < +00 for all t > 0. 

• ^Viep -G*p)£ L2([0, T] X M2) for all T > 

such that the following energy identity is satisfied 



E[p{t)]+ [ I p\V{ep-G*p)\'dxdt = E[pQ]. 

Jo JR'* 

In particular, the equation (2.1) produces the following regularizing effect. 



(2.3) 



Lemma 2.2 (Regularity of H P steady states). Let pQ G L^ nV. Then, the corre- 



sponding solution p{t) to (2.1) satisfies 

j p |Vpp dx < +00, 



(2.4) 



for almost every t > 0. In particular, let p be an L^ CiV steady state to (2.1), then p 
satisfies (2.4) and p ^ G"^ on supp[;9]. 



Proof. Due to (2.3), the quantity 



e j p\Vp\^dx-2£ j pVp-VG* pdx + J p\VG*p\'^dx 

is finite for almost every t > 0, and therefore, in view of Cauchy-Schwarz inequality, we 
have 

ly piV/jpda; - C(e) J p\VG*p\'^dx + J p | VG * < +oo. 
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and thanks to the smoothness assumptions on G we have the assertion (2.4). Let p be 



a steady state, then p satisfies (2.4) too. This implies in particular that Vp is almost 
everywhere finite on M^. The energy identity (2.3) implies then 



p\V{ep-G* p)\^ = 
for almost every j; G R^. This means that 

ep — G * p = constant 

almost everywhere on every connected component of the support of p. By convolution 
with standard mollifiers, one can easily see that ep — G * p = C for a given C depending 
on the connected component of supp[/9]. Since G is C^, this easily implies p G on 
supp[p]. □ 



Corollary 2.3 {Id regularity). Let p be an L^H'P solution to (2.2) in one space dimen- 
sion. Then p is continuous on M. 



Proof. Apply the result in Lemma 2.2 to the case d = I. Since £p^/^ G L^, the one 



dimensional Sobolev embedding implies that p^/^ is continuous. □ 

3. Stationary solutions in multiple dimensions 

In this section we state the necessary and sufficient conditions on e and such 
that there exists non trivial steady states 

pV{ep-G*p) = (3.1) 

in the set Pi V. During our work, we realized that J.Bedrossian has obtained similar 
results in [2j , based on ideas and strategies developed in [38] and [ID] . In order to simplify 
the coverage of the paper, we shall state the result in [21 Theorem 1] and prove all other 
results. Notice that the critical case £ = ||G||j;^i was not covered in [2]. 



Let us start by focusing on the interplay between the solutions to (3.1) and the varia- 
tional calculus on the energy functional 

1 



E[p] :=£ p {x)dx - - / G{x - y)p{y)p{x)dydx. 
In the next proposition we prove that being a minimum for the entropy functional is a 



sufficient condition for being a solution to (3.1) 



Proposition 3.1 (Stationary solutions via energy minimization). Let p G L^(M'^) be a 
minimizer for the energy functional 



E[p] ■■=\ f p{ep-G*p)dx 



on V . Then 

pV (ep-G* p) = a. e. in M'^. 

Proof. Let V G C^(M'^) be an arbitrary vector field and let u{x,s) be a local solution to 
the continuity equation 

dsu{x, s) + div(n(2;, .s)V{x)) = 

with initial datum 

u{x, 0) = p{x) 
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with p being the minimizer for E given in the hypothesis. Such a u can be constructed 
by solving the characteristic ODE 

j-X{x,s) = V{X{x,s)) 

coupled with the initial datum 

X(x,0) = X 

locally in s = 0, with the local solution X{x,s) being C^, and by taking u{x,s) := 
s))(|/o](x, s), i. e. s) defined via 

x)u{x,s)dx = j (p{X{x,s))p{x)dx, for ah (p £ Cl{R'^) 

(cf. for instance Chapter 8, Lemma 5.5.3]). For all s in the interval of existence of u 
we have 

J u{x,s)dx = l, u{x,s)>0 a.e. 

and therefore the map s i— t- E[u{-, s)] has a local minimum at s = 0. Hence 

d 
ds 

= — J (eu — G *u) div{uV)dx\s=o = J P^i^P — G * p) ■ Vdx 
and replacing V with —V we obtain 

0> J pV{ep -G*p)- Vdx 

and therefore 

j pV{ep -G* p)- Vdx = 0, for an arbitrary V G GI{W^) 

which is the desired assertion. □ 
Let us now compute the first and the second order Gateaux derivatives of E. 



< ^E[u{-, s)]|s=o = I (eu-G* u)dsudx\s=o 



Lemma 3.2. Let p £ L nV be a solution to (3.1 ). Then, p is a stationary point for the 
energy functional E. Moreover, the second order Gateaux derivative of E on p satisfies 

d^ f f 

-tzkEIp + 6v]\s=o = e / v'^(x)dx— / v(x)G * v(x)dx, (3.2) 
dd"^ ./rod ./ 



for all V = diy{pV) and V G G^iJ 
Proof. 



Suppose p £ L nV satisfies (3.1). Let us compute 

hm-(E[p + 5v]-E[p]) 
s^o 
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with V = div{pV) for an arbitrary vector field V £ C^, which imphes J v{x)dx = 0. We 
obtain 

\e[p + 5v]-E[p]) 



6 

{p + ^vfdx - ^ I P^dx 

■^^ J supp{p+Sv) Jsupp(p) 

- I {p + 6v)G * {p + 6v)dx + I pG*pdx. 

Therefore we easily get 

lim \ (E[p + 6v] - E[p]) = [ v(ep -G* p)dx 
<5->-0 J 

= j <ii\{pV){£p -G* p)dx = - j pV ■ V{ep -G* p)dx. 

Therefore, p\s a. stationary point for E under the constraint J pdx = 1. The computation 
of the second derivative of -E on p yields 

^E[p + dv] = -7^2 7^ {p + 6vfdx 



dS^ 

d^l 

d6^2 



supp(p+(5ii) 

£ / v"^dx — / vG * vdx 



' supp(p+(5t)) 

{p + 5v)G * {p + Sv)dx 



which is independent on 6 and therefore it is valid also on 5 = 0. □ 

Before we start analyzing the existence or nonexistence of steady states, we introduce 
a very simple technical lemma which will be very useful in the sequel. 

Lemma 3.3. Suppose p £ L"^ nV is a solution to (3.1) having connected support. Then 



ep{x) = G{x-y)p{y)dy + G 

Js\ipp[p\ 

for all X G supp[p] with G = 2E[p]. Moreover, in case supp[/?] has infinite measure, then 
C = E[p\ = 0. 

Proof. It is immediate from (|3.1|) that 



ep{x) = / G{x-y)p{y)dy + G (3.3) 

Jsupp[p] 



for all X G supp[yo] for a certain constant G. Then, we multiply (3.3) by p{x) and integrate 
over supp[/9] to obtain 

p^{x)dx= / / G{x — y)p{y)p{x)dydx + G, 

supp[p] J supp[p] J supp[p] 

where we have used that p has unit mass. It is therefore clear that G = 2E[p]. Suppose 
now that supp[/9] has infinite measure. Suppose by contradiction that C / 0. Let 
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{xk} C supp[/9] be a sequence of points such that \xk\ — )• +00. We have, for ah k, 

£p{xk) - I G{xk - y)p{y)dy = C 

Jsupp[p] 

and therefore the same expression should hold in the limit k — >• +00. Now, the assump- 
tions on G imply that the integral 

/ G{xk - y)p{y)dy 

Jsupp[p] 

converges to zero as A; — )• +00. This is due to Lebesgue's dominated convergence Theorem. 
Therefore, the term p{xk) has a limit C as A; — t- +00. Such limit is the same for all 
diverging sequences of points {x^} C supp[/9], which means 

lim p(x) = C. 

a;£supp[p],|a;|— >+oo 

Now, since supp[/9] has infinite measure, then C 7^ implies that p is not integrable, 
which is a contradiction. Therefore C = 0. □ 



3.1. Non existence of nontrivial steady states for e > 1. We start by covering the 
case £ > \. Here, there exist no nontrivial steady states, as it follows from the following 
simple lemma. 



Lemma 3.4. Let e > 1. Then, there exists no stationary solutions to (3.1) in the space 

L'^nv. 

Proof. We first prove that there exists no minimizer for E[p] under the mass constraint 
j p = 1 and p >0. To see this, we use Young inequality for convolutions as follows 

E[P] = 1 1 P'd^-lj pG*pdx>'- I p^dx-&^ I p'dx = '-^ I p'dx (3.4) 

with e — 1 > 0. Moreover, we have the simple estimate E[p] < C||p||^2- Take a family 
of functions px{x) > such that / px{x)dx = 1 and / p\{x)dx — )■ as A — )• +00. To 
construct such a family, we just take a fixed L^(M'^) function p ^ and rescale it by 
p\{x) = \~'^ p{\~^x). For such a family we therefore have 

E[p\] — )• 0, as A — )• 00. 

Therefore, it is impossible to have a minimizer p^ for E[p\ in the set {p G : = 1} 



because (3.4) would imply that E[poo\ > and we would necessarily have < E[px] < 
E[poo] for A large enough. 

Now we prove that there exist no steady states. Suppose by contradiction that p is 



a steady state. Then, due to Lemma 3.2 /) is a stationary point for E. Moreover, the 
formula (3.2) implies that the functional E is strongly convex, and therefore admits only 
one stationary point, which coincides with its global minimizer. But this contradicts the 
non existence of a global minimizer proven above. □ 
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3.2. The critical case e = 1. We aim to solve 

= div(pV(p-G*p)). (3.5) 

We shall prove that no L? r\V steady states exist in this case. 

Theorem 3.5 (Non-existence of nontrivial steady states for e = 1). There exists no 
solutions to (3.5) inL'^CiV. 

Proof. From Cauchy-Schwartz inequality we know that 

j^^pG* pdx < ||p||L2(Rd)||G*p||^2(B;d) 

and the equality in the above formula holds if and only if p and G * p are proportional. 
In terms of the functional E this means that 

E[p] > for all p G L'^{R'^) n P(M'^). 

3.4, we have the estimate E[p] < C||p||^2) and using once again the family 



As in Lemma 
p\ of Lemma 



3.4 we see that 'mfp^i2 E[p] = 0. Assume by contradiction that there 
oo . Then, due to the result in Lemma 13.21 and in view of 



exists a stationary solution p, 
Cauchy-Schwartz inequality, the second order derivative of E is nonnegative everywhere. 
Hence, the functional E is convex and therefore poo is a global minimizer for E under the 
constraint p ^ CiV. Then, we must have E[poo] = 0, which means that p and G * p 
are proportional, i. e. there exists a constant A G M+ such that 

Pooix) = XG * poo{x) (3.6) 

almost everywhere on M'^. Integrating (3.6) over M*^ yields 

= A 



1 



M\Gh 



and hence 

almost everywhere on 



Poo{x) = G* poo{x) 
l'^. We can then apply the Fourier transform 



(3.7) 



fiO 



-2-«-€/(x)(ix 



to both members of the equation ( |3.7| to obtain 
We have 

|G(OI < / \Gix)\dx = l. 

Moreover, since G is even, then G{(,) < 1 for all 7^ 0. In order to see that, write 

d „ d 



i s\.n{2TiXk^jS))G{x)dx , 



G{0= [ ne-2™^«'=G(x)(ix= / rr(cos(27rx,,efc) 

then G being even easily implies that only real valued contributions survive in the above 
integral; such real valued contributions are of the form 

fh,k{x,^)G{x)dx 
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where the functions fh^k are such that \ fh,k{x,S,)\ < 1 and \ fh,kix,S,)\ < 1 for j; ranging 
on a set of positive measure. Therefore, we have proven that 

p^{0 = for all 

and p^{0) = 1. This implies that p{x) = almost everywhere, which contradicts the fact 
that p has unit mass. □ 

3.3. Stationary solutions for e < 1. Let us now provide a minimizer for the entropy 
functional in the case e < 1, which implies the existence of a nontrivial L'^ CiV steady 



state for (3.1 ) in view of Proposition 3.1 Such result is proven rigorously in [2J Theorem 



1], which we recall here. 

Theorem 3.6 (Existence of minimizers, [2J). Let e < 1. Then, there exists a radi- 
ally symmetric non-increasing minimizer p ^ V H L^(M'^) for the entropy functional E 
restricted to V with p ^ 0. 

We send the reader to [2j for the details of the proof, which is based on a sort of sub- 
additivity property needed to provide suitable compactness of the minimizing sequence, 
cf. |2t Lemma 2]. For the sake of clarity, we shall still provide the simple proof of the 
fact that global minima of E under mass constraint are strictly negative, which forces 
the minimizer to be non zero. 

Lemma 3.7. Let e < 1. Then, infpg-pp|^2(]g<i) E[p\ < 0. 
Proof. We consider the family ax £ L^ CiV 



For e < 1 we have 



^[^^^ = ^ - sl^ /, /, - y)^^^^ 

£ 1 /"^ . . 1 ^ '■^ 



4A 4A./_,'^(^><''- = 4a(,^-I,«<^'''^ 



and since 



.A 

J G{z)dz — )• 1 as A ^ +oo, 
we easily obtain that there exists a A such that E[a-^ < 0. □ 

4. Stationary solutions in the l-d case 

In this section we prove the main result of our paper, namely that nontrivial stationary 
solutions (which always exist in the case e < 1) in one space dimension with fixed mass 
and center of mass are unique. We shall first provide certain necessary conditions on the 
steady states and then prove that they are unique under such conditions. The main tool 



in this procedure is the use of the strong version of Krein-Rutman Theorem 4.10 

We start with a necessary condition on the steady states which deals with a property 
of their support. 
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Lemma 4.1 (Steady states have connected support). Let p be a stationary solution to 
(3.1) in one space dimension, namely 

pdx{ep - G * p) = a.e. onR. (4.1) 

Then, supp(p) is a connected set. 



V{x) :-- 



Proof. Let p solve (4.1). Let us first assume that p is compactly supported. Suppose that 
supp(/9) is not connected. Accordingly, let [a, b] be a non trivial interval such that 

p{x) ^0 if X < a 

p{x) = ifa<x<b (4.2) 

p{x) ^0 if X > b. 

Let us introduce the velocity field 

— 1 if X G (— oo, a) 
1 if X G (6, +oo) 

and let V € C^(M). Let u{x,s) be a local solution to the Cauchy problem for the 
continuity equation 

(d,u + d,{uV) = o 

\u{x,0) = p{x). 

Let us compute the evolution of the energy E along u at the time s = 0: 
d_ 
ds 

Then, by definition of V we have 

£ / pvd^p = - / d^p^dx + - / vd^p^dx - - / d^p^dx = (4.3) 

Jr ^ J-oo Ja Jb 

because = on [a, b] and p = 0onx = a,6 and at ±00. Hence, we have 

/4-00 pa r+oo 

pVdccG *pdx = - / pG' *pdx+ / pG' * pdx. (4.4) 
00 J — 00 J b 

We compute 

/a pa ra r+oo 

/ p{x)G'{x-y)p{y)dydx+ j / p{x)G' {x - y)p{y)dydx, 

the first term on the above right-hand side is zero since G' is odd and the integration 
domain is symmetric in x and y. Since ^'(z) > as 2: < 0, we have for the second term 

p{x)G' {x — y)p{y) > on (x, y) G (—00, a) x (6, +00). 

In a similar way one can prove that 

+00 r+oo pa 

pG' * pdx = / dx dy p{x)G' {x — y) p{y) 



-E[u{s)\\s=o = j Us{£u{x, s) — G * u{x, .s))dx\s=Q = j pVdx{£p — G * p)dx = 0. 



with the integrand p{x)G'{x — y)p{y) < on the integration domain. Therefore, (4.4) 
implies that 

p{x)p{y) = on {x < a} n {y > b}. (4.5) 
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We have thus proven that, whenever (4.2) holds, then (4.5) has to be necessarily satisfied. 
Let A, B be two nonempty connected components of supp(/9) and let [a, /?] be the maximal 
interval such that 

a < 6, for aU a G A, h & B 
p = on [a, 13] 
a > a, for all a G A 
f3 <b, for aU b £ B. 
Then, p{x)p{y) = for all such that 

X < a, y > /3, 

which implies that either A or B cannot be in the support of p, and that is a contradic- 
tion. In order to generalize the proof to a stationary solution p which is not compactly 
supported, one can cutoff p to have compact support in such a way that the norm of 
the compactly supported approximation is arbitrarily close to the norm of p. Then, 
the estimate 

e j \pdxp\dx < ||G||ii||/9||^2 

implies that the integrals in (4.3) converge at infinity, therefore all the above computations 
are valid up to an arbitrary difference which vanishes in the limit. □ 

Remark 4.2. In the case supp(G) = [— ii, B\ one can use the same strategy as in Lemma 



4.1 to prove that, given two connected components A, B of supp(/9) one has dist(^, B) > 



2R. The proof is a straightforward generalization of the above arguments, and it is 
therefore left to the reader. 

We now exploit a standard symmetric rearrangement technique to prove that the min- 
imizers of the energy are symmetric and monotonically decreasing on a; > under the 
constraint of zero center of mass, cf. |38t I39j. 

Proposition 4.3. Let p^o be a minimizer for the energy 

E[p] = 7; I P^{x)dx- \ f 1 G{x - y)p{x)p{y)dydx 



2 7 2 „ „ 

under the constraint that the center of mass is zero. Then, poo is symmetric and mono- 
tonically decreasing on x > 0. 

Proof. We have to prove that the energy decreases strictly when a function u which is 
not symmetric and decreasing on x > is replaced by its symmetric rearrangement 

u*{x) = sup{t > : meas({n > t}) > 2 |x|} . (4.6) 

For every exponent p > 1 the following holds: 

[ {u*)Pdx = [ {uYdx, (4.7) 

therefore the part of the energy is invariant when passing from u to u*. As for the 
interaction energy, we recall the well known Riesz's rearrangement inequality, see e. g. 



/ / 

JR'' JR" 



f{x)g{x-y)h{y)dydx< / / f* {x)g* {x - y)h* {y)dydx, (4.8) 
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which holds for all nonnegative functions /, g, h vanishing at infinity. Moreover, if g is 
strictly decreasing on a; > and symmetric, then equality in (4.8) holds if and only if 
/(x) = f*{x — xq) and h{x) = h*{x — xq) for some common xq. Apply such a theorem to 
our case, using G* = G and the fact that u is not symmetric up to translations. □ 

Let us rephrase Lemma 3.3 in the one-dimensional case. 
Lemma 4.4. Let p be a Id steady state, i. e. 

ep = G * p + C on supp[/9] 
for some C G M. Then, C = 2E[p]. 



Proof. The support of p is connected in view of Lemma 4.1 , therefore Lemma 3.3 applies. 

□ 

Lemma 4.5. Let p £ V Ci L"^ and let xq G M. Let px^ be defined by 

Pxoix) ■■= p{x + Xq). 

Then, E[px,;\ = E[p]. 

Proof. It follows by direct computation of the energy and by change of variable under 
the integral sign. □ 

Lemma 4.6. Let p be a steady state with e < 1. Then, the support of p is compact. 



Proof. We know from Lemma 4.1 that the support of /) is a connected set. Suppose that 
supp(/9) is not bounded. That means that supp(/9) is of the form (—00,6) (6 possibly 



+00) or (a, +00) (a possibly (—00). Assume first supp(p) = (a, +00). Then, Lemma 4.4 
implies 



r+00 

2E[p] = ep{x) - G{x- y)p{y)dy = 

J a 



for all x G (a, +cxd). Now, there are two possibilities: either a = — oo or o > — oo. In the 
latter case, evaluation on a; = a implies 

r+co 

= ep{a) = G{a- y)p{y)dy 

J a 

which is a contradiction because the integral on the right hand side is strictly positive 
in view of supp(G) = M. In the former case a = — oo we have then supp(/9) = M, which 
implies 

r+oo 

ep{x) = I G{x - y)p{y)dy 



for all a; G M. We can therefore integrate over M to obtain 

e= IIGIlii = 1 

which is a contradiction. The same proof can be produced in the case supp(/9) = (— oo, b). 
Therefore, the support of p can only be a bounded interval. □ 

Lemma 4.7. Let p be a steady state. Then there exists a symmetric steady state p such 
that 

E[p\=E[p]. 
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Proof. From Lemma 4.1 and Lemma 4.6 we know that supp[/o] = (a, b) for some a, 6 G M. 
For a given x G (a, 6) we have 

£p(x) = G* p{x) + C (4.9) 
for some C G M. Evaluation on a; = a and x = b gives 

C = - f G{a- y)p{y)dy = - f G{b - y)p{y)dy. 

J a J a 

Let p{x) = p{x + xq) with xq = {a + b)/2. Then p is still a steady state and it satisfies 
E[p] = E[p] thanks to Lemma 4.5 Moreover, the support of p is symmetric. Let us 
introduce 

Kx) ■■= ^{p{x) + p{-x)). 
Clearly, supp[p] = supp[p] and we have, for all x G supp[/)|, 

^Pix) = ^{p{x) + pi-x)) = |(p(x + xo) + pi-x + xo)) 

1 1 

= 2 + xo- y)p{y)dy + - G{-x + xq - y)p{y)dy + C 

= - G{x- z)p{z)dy + - / G{-x - z)p{z)dy + C 

^ J{a~b)/2 ^ J{a~b)/2 

1 fip-a)/2 ^ r{b--a)/2 

G{x - z)p{z)dy + - / G{x- z)p{-z)dy + C 

(a-b)/2 ^ J{a-b)/2 

{b-a)/2 I r{b~a)/2 

G{x - z)-{p{z) + p{-z))dz + G = / G{x - z)p{z)dz + C 

(a-b)/2 2 J(a-b)/2 

where we have used the symmetry of G. The above computation shows that p has the 



same energy as p in view of the results in Lemma 4.5 and Lemma 4.4, □ 



Lemma 4.8 (Support of a minimizer). Let poo be a global minimizer to E. Let p be a 

steady state such that 

meas(supp[poo]) < meas(supp [/?]). 
Then p is also a minimizer. 

Proof. Assume first that we are in the special case supp[/3oo] ^ supp[/9]. Let us compute 
the second variation of E around the minimizGr ciloiig ttL6 direction — p. 

-^E[poo + S{p - poo)]\5=0 



e y (p - Poofdx - J J G{x - y){p{x) - Poo{x)){p{y) - poo{y))dydx 
2E[p\ + 2E[p^] -2 f p^{ep-G* p)dx 

J sunn f 0^1 



/supp[p, 

2E[p\ + 2E[pJ^ - ^E[p] 
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where the last step is justified by the fact that supp[/3] C supp[/>oo]- Therefore, since p^o 
is a minimizer, the second derivative above is nonnegative, i. e. 

< ^^[Poo + - Poo)]\s=o = 2{E[p^] - E[p]), 

which yields E[p\ < E[poo]. Since poo is a minimizer, then so is p. In the general case 
in which supp[/Ooo] ^ 

supp[/3], consider a translation Pxq{x ) = p{x - xq) in such a way 
that the support of p contains the support of poo- Since the energy is invariant after 



translation in view of Lemma 4.5, the assertion is proven. □ 



We are now getting closer to the proof of our uniqueness result. Let us recall the 
following important theorems, see e. g. [23j and the references therein. 

Theorem 4.9 (Krein-Rutman Theorem). Let X be a Banach space, let K C X be a 
total cone, i. e. such that XK C K for all X> and such that the set {u — v, u,v £ K} 
is dense in X . Let T be a compact linear operator such that T{K) C K with positive 
spectral radius r(T). Then r{T) is an eigenvalue for T with an eigenvector S i^T \ {0}. 

An important consequence [23j of the Krein-Rutman theorem, which will be extremely 
useful in the sequel, is the following 

Theorem 4.10 (Krein-Rutman Theorem, strong version). Let X be a Banach space, 
K C X a solid cone, i. e. such that XK C K for all X > and such that K has a 
nonempty interior Kq. Let T be a compact linear operator which is strongly positive with 
respect to K, i. e. such that T\vl\ G Kq if u £ K. Then, 

(i) The spectral radius r[T) is strictly positive and r{T) is a simple eigenvalue with an 
eigenvector V € Kq. There is no other eigenvalue with a corresponding eigenvector 
V £ K. 

(a) \X\ < r[T) for all other eigenvalues X ^ f'{T). 

We shall now prove the uniqueness of symmetric steady states with unit mass which 
are monotonically decreasing on the positive semi-axis in the case e < 1. We already 
know that under the above assumptions we can write, for x £ [—L,L] = supp[/j]. 



ep{x) 



j'^^G{x-y)p{y)dy + C, C = 2E[p]. (4.10) 



Taking the derivative w.r.t x £ [—L, L] we obtain 

, d d / 

£p\x) = — J ^ G{x - y)p{y)dy = —G * p{x) = J ^ G{x - y)p{y)dy. 

The symmetry of p and G implies, for x £ [0, L], 

pL pL pL 

ep\x) = - G{x+y)p'{y)dy+ G{x-y)p' {y)dy = [G{x - y) - G{x + y)] p' {y)dy. 
Jo Jo Jo 
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Assuming that p G C^([— finding a steady state with the above assumptions is 
equivalent to find p on [0, L] such that 

p{L) = 0, 

— p{x) = u{x), xG[0,L], 

u > 0, and u solves eu = H{x, y)u{y)dy, 

Jo 

H{x,y) = G{x-y)-G{x + y). 
To convince ourselves about that, integrate 

-ep'ix) = - [\g{x -y)- G{x + y))p'{y)dy 



over the interval L] for some ^ G [0, L). Then p{L) = and integration by parts imply 

-ep{i) = - I dx [ {G{x-y)-G{x + y))p'{y)dy 
Jo 

dx[{G{x -L)- G{x + L))p{L) - (G(x) - G{x))p{0)] 

r-L 



+ [ dx [ {-G'{x-y)-G\x + y))p{y)dy 
Js, Jo 



L 







p{y)dy / {-G'{x - y) - G'{x + y))dx 



: r p{y)[-G{L -y)- G{L + y) + G{i - y) + G{i + y)] 
Jo 



which implies, by the symmetry of G, 
ep{x 



[ G{x - y)p{y)dy + G, C = - [ {G{L - y) + G{L + y))p{y)dy. 

J-L Jo 
For further reference, we introduce the operator 

Gl[p\{x) := [G{x -y) + G{x + y) - G{L - y) - G{L + y)] p{y)dy (4.11) 
Jo 

on the Banach space 

Yl:={p^C{[Q,L]) : p{L) = 0} . 
In order to simplify the notation, we also define the following operator 

Ul[u]{x) := I H{x,y)u{y)dy = {G{x - y) - G{x + y))u{y)dy. 
Jo Jo 

Proposition 4.11. For a fixed L > there exists a unique symmetric function p G 
C^([— L,L]) with unit mass and with p'{x) < on x > such that p solves (4.10) for 
some £ = £{L) > 0. Such function p also satisfies p"{0) < 0. Moreover, £{L) is the 
largest eigenvalue of the compact operator Ql on the space Banach Yl and any other 
eigenfunction ofQi on Yi with unit mass has the corresponding eigenvalue e' satisfying 



QUADRATIC DIFFUSION EQUATIONS WITH LONG-RANGE ATTRACTION 17 

Proof. Since G is decreasing on the half-line [0, +C!o) we get 

H{x, y) = G{x — y) — G{x + y) > 0, on x,y > 0. 
Consider now the Banach space 

XL = {fGGH[0,L]): /(0) = 0} 

endowed with the norm 

ll/IUi = II/IIl°=([o,l]) + II/'IIl°°([o,l])- 

It can be easily seen that the set 

K:={feX: />0} 

is a solid cone in X. Indeed, any function f G K with /'(O) > is in the interior of K. 
Moreover, for a given u G K, we have 

niiuKx) = [ Hix,y)u{y)dy>0 
Jo 

for all x S [0,L] and 

^lN(0)= ["^ H{0,y)u{y)dy= [\G{-y) - G{y))u{y)dy = 0. 
Jo Jo 

Therefore Ti is a positive operator in the sense provided by the definition of the cone K. 

Indeed, we can prove that V. is strongly positive, i. e. for a given u £ K, 7i[u] belongs to 

the interior of K. In order to see that, for a u £ K \ {0} compute 

CHlM)'(O) = f\G'{-y) - G\y))u{y)dy = -2 G\y)u{y)dy > 0, 
Jo Jo 

and therefore T-Ll[v\ belongs to the interior of K. Hence, we can apply the stronger version 
of Krein-Rutman theorem 4.10, which implies the existence of a simple eigenvalue e > 
equal to the spectral radius oiT-Li. More precisely, there exists a family of solutions u to 

eu = 'Hl [u] 

generated by one given nontrivial element u in the interior of K. This implies that the 
corresponding set of symmetric and monotone p solving (4.10) satisfies 

fL f'L f'L 

p{x) = p{x) - p{L) = - p{y)dy = / u{y)dy = a u{y)dy, 

J X J X J X 

for a > 0. We choose a as follows 



and we obtain that p has unit mass on [— L, L]. It is clear that p'{x) < for a; > 0, 
p'(0) = 0, and p"{0) < 0. In view of the statement (i) of theorem 4.10, there exists no 
other eigenvalues to Hl with eigenvectors in K besides the one e with eigenfunction u, 
and all other eigenvalues e' with eigenfunctions in Xl satisfy |e'| < e. □ 

The eigenvalue e (which coincides with the spectral radius oi T-Ll) can be considered 
as a function of L, namely e = £{L). The behavior of such function is established in the 
next proposition. 
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Proposition 4.12 (Behavior of the function £{L)). The simple eigenvalue e{L) found in 
Proposition \4 - 1 i is uniquely determined as a function of L with the following properties 

(i) £{L) is strictly increasing with respect to L 

(ii) HmL^+oo £{L) = 1 

(iii) £{0) = 0. 

Proof. In order to prove (i), let us consider the equation 

e{L)uL{x) = Hl[ul]{x) = H{x,y)uL{y)dy, xe[0,L], 

Jo 

where ml is the unique eigenfunction obtained in Proposition 4.11 We multiply the above 
equation by ul{x) and integrate over [0,L] to obtain 

e(L) / UL{x)'^dx= I 'HL[uL\{x)u[x)dx. 
Jo Jo 

Recall that the eigenvalue ul satisfies ^^(O) = and, for x G (0, L], 

'^l{x) = 1 H{x,y)u{y)dy > 



e{L) Jo 

since H{x, y) = G{x — y) — G{x + y) > for all y G [0, L] under the assumption x > in 
view of the strict decreasing monotonicity of G on a; > 0. For a general L £ (0, +oo) and 
a 5 > (small enough) we have 

rL+S 



h := e{L + 6) / u'j^^^{x)dx — £{L) / u'j;^{x)dx 
Jo Jo 



L+S i-L 

T-l'L+5[uL+5]{x)uL+5{x)dx - 'HL[uL\{x)uL{x)dx =: h- (4.12) 
'0 Jo 

We analyze the two terms Ii and I2 separately. Ii can be expanded as follows: 

r-L+5 

.2 



/i = {e{L + 5)- e{L)) [ 

Jo 

+ e(L) / {ul+s{x) - ul{x)){ul+5{x) + UL{x))dx + e{L) u\{x)dx. 
Jo Jl 



UL+5{x)dx 

fL+6 i-L+S 



I2 is given by 



h = I {T-Ll+5[uL+5]{x) -HL[uL]{x))uL+delta{x)dx 





+ I 'HL[uL]{x)uL{x)dx + I niiuLjix) (uL+six) - ul{x)) dx 

L Jo 

L+5 

{UL+i [uL+s] {x) - Ul [ul] [x] ) ui+deita {x)dx 



f-L+S f-L+S 

ul{x) (uL+six) - ul{x)) dx 



+ e{L) / u'l{x)dx + £{L) 

Jl Jo 
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and on substituting Ii and I2 in (4.12) we can cancel some terms and obtain 

{e{L + 6) - e{L)) u\j^^{x)dx + e{L) l {ul+s{x) - UL{x))uL+s{x)dx 

Jo Jo 

f-L+S 

{V.L+s[uL+5]{x) - V-lIulKx)) UL+5{x)dx 



fL+<5 rL+S 

H{x, y)uL{6){y)uL+5{x)dydx 

10 JL 

+ / / H{x,y){uL+s{y) -UL{y))uL+5{x)dydx 
Jo Jo 

r-L+S f-L+S 

H{x, y)uL{y)uL+5{x)dydx 

lo JL 

f-L+S 

UL+s{x){uL+six) - UL{x))dx 



+ e{L + 5) [ 
Jo 



where we have used the definition oH-Ll and the property H{x, y) = H{y, x). By suitably 
expanding the term on the left hand side in the above identity, we obtain 

i-L+S i-L+S j-L+S 

{e{L + 5) - e{L)) / ULix)uL+s{x)dx = / / H{x,y)uL{y)uL+s{x)dydx 

Jo Jo JL 

t-L+S 

= e{L + 6)J^ UL+s{y)uL{y)dy 

and the positivity property of ul implies that 

eiL + 6) > e{L), 

which proves (i). 

Let us now prove (ii). Assume by contradiction that 

lim e{L) = eo < 1- 



Let e G {eo, 1). We know from Theorem 3.6 that there exists a minimizer for the 
energy E with zero center of mass. We also know that the support of ps is compact 
from Lemma 4.6 From Proposition 4.3 we know that ps is symmetric and monotonically 
decreasing on x > 0. Therefore, ps is the unique eigenfunction with unit mass provided 
by Proposition 4.11 and the support of pe is [— L,L] for some L > 0. Therefore, the 
corresponding eigenvalue should be e{L) < eq, which is a contradiction since e and Eq are 
two different eigenvalues with the same eigenfunction. 

Let us prove (iii). By letting L\0 one has that the operator Tii'^s the zero operator, 
and therefore e(0) should be the eigenvalue of the zero operator, which can only be zero. 

□ 

We are now ready to prove the main result of this paper. 
Theorem 4.13. Let e <1. Then, there exists a unique p £ solution to 

pdx{ep- G*p) = 0, 
with unit mass and zero center of mass. Moreover, 

• p is symmetric and monotonically decreasing on x > 0, 
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p G C2(supp[p]), 

supp[p] is a hounded interval in M, 

p has a global maximum at x = and p"{0) < 0, 

p is the global minimizer of the energy E[p] = | J p'^dx — | / pG * pdx. 



Proof. We know from Theorem 3.6 that there exists a minimizer poo with miit mass 
and zero center of mass, which is symmetric and monotonicaUy decreasing on x > in 
view of Proposition 4.3 and compactly supported on a certain [—L, L] in view of Lemma 



4.6 From the results in Propositions 4.11 and 4.12, we know that there exists a unique 
steady state with such properties, because the correspondence e = e{L) is one-to-one. 
So, the only possibility to violate uniqueness of steady states with unit mass and zero 
center of mass is to have a steady state which violates either the monotonicity property 
or the symmetry. Suppose first that there exists a steady state with zero center of mass 
p which is not symmetric, it is not restrictive to assume the support of p is [—L',L'\. 



Then, we know from Lemma 4.7 that it is possible to construct a symmetric steady 
state p with the same energy of p and with the same support of p. Now, there are two 
possibilities: either p is a minimizer or not. In the former case p is also a minimizer and 
this is a contradiction (a minimizer is symmetric). In the latter case, the support of p is 
strictly contained in the support of poo in view of Lemma 4.8 and p is not monotonicaUy 
decreasing on a; > because otherwise it would be the unique minimizer provided before. 
Therefore, with the notation of Proposition 4.11 —pi' is an eigenfunction for %ii in the 



space which is not belonging to the solid cone K. Therefore, the stronger version of 
Krein-Rutman's Theorem 4.10 and the fact that e{L) is increasing imply that L' > L, 



since p is an eigenfunction outside the solid cone -fC, and it therefore should have an 
eigenvalue strictly less than e{L'). This implies that e{L) < £{L') and therefore L < L' . 
Now this is a clearly a contradiction because we said before that the support or p is 
strictly contained in the support of poo, so L > L' . The case in which p is symmetric but 
not monotonic on x > can be covered by repeating the same argument above (assume 
p = p!). □ 

Corollary 4.14 (Concavity of p for small e). There exists a value Eq € (0, 1) such that, for 
all e € (0,eo) ^/^e corresponding stationary solution provided in Theorem 4-13 is concave 
on the whole interval [0,L]. 

Proof. We can differentiate twice w.r.t x in 

ep{x) = / G{x - y)p{y)dy + C 



to obtain 



ep"{x) 



G"{x - y)p{y)dy 



for all X G [— L,L]. Therefore, G" is evaluated on the interval [— 2L, 2L] in the above 
integral. We know from Proposition |4.12 that L is a monotonicaUy increasing function 
of e with lim^x^o L{e) = 0. Since G"(0) < 0, and G G C^, then there exists Lq > such 
that G" < on [— 2Lo,2Lo]. Let eo be the eigenvalue in K corresponding io L = Lq. 
Then, the eigenfunction p is concave on its support. □ 
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Remark 4.15 (The case meas(supp(G)) < +00). It is also interesting to consider the case 
with the support of G being bounded, i.e. supp(G) = [—g,g\ with G being symmetric and 
monotone on [—g, 0]. Then most techniques of the paper remain vahd, in particular when 
considering simply connected stationary states, whose existence can be shown along the 
lines of the above arguments. With the same arguments as in the proof of Lemma 4.1 one 
can show that the distance between two connected components of a stationary solution 
is at least 2g. This means that stationary solutions consist of a countable number of 
connected components, which are not influenced by any other connected component (the 
distance is larger than the kernel range). Thus, each connected component is a stationary 
solution by itself. Thus, in particular we can prove the existence of infinitely many 
stationary states. The remaining step would be to characterize uniquely the behavior 
of each connected component as an energy minimizer. All arguments above apply to 
the case of simply connected solutions in the case of a finite range kernel G, except the 
application of the strong version of the Krein-Rutman theorem. If we cannot guarantee 
that the support of the global minimizer is contained in [— f) §]' there can be elements in 
the nonnegative cone such that the convolution with G is not strictly positive. Although 
we strongly believe that such a smallness of the support holds at least for e small, we so 
far did not succeed in proving such a result. 

5. Numerical Results 



In the following we show numerical simulations for the evolution equation (2.1). We 
discretize the equation using an explicit Euler scheme and finite difference methods. In 
one dimension, we partitionate the domain VL = [a, h] using an equidistant grid with n + 1 
grid points a = xq < xi < ... < Xn = h and step size h = [b — a)/{n + 1). Furthermore 
we use the following finite difference scheme: 



dt --vr^--^^-r^ - r« 

with forward and backward difference quotients 

A+i - Pi , _ Pi ~ Pi-^ 

JJ+Pi- ^ , U^pi- ^ . 

The time step size dt has to be chosen appropriately, to guarantee stability. In a first 
example we consider an interaction potential G{x) = — \= exp {—\ (^^^)^ ) with mean 



/2-K 

jjL = Q and variance o" = 1, which fulfill the conditions (l)-(8), HGH^^i = 1 and models a 
wide range attraction. For this kernel we present in Figure [T] the solutions for the station- 



ary equation (2.2), which means we calculated the largest eigenvalues and corresponding 



eigenvectors of the operator Ql defined in (4.11) for different L. The largest eigenvalues 
e = e(-L) are presented in Figure [T] (a). Like mentioned in Proposition 4.12 e(L) is strictly 
increasing with respect to L and furthermore lim^^-^oo ^{L) = 1. The corresponding eigen- 
functions with unit mass are presented in Figure [T] (b) and (c). We proved in Corollary 



4.14| the concavity of p for small e. To better clarify the situation, we illustrate in Figure 
1] (b) the eigenfunctions for L S (0,1], that means for e < 0.05. For a certain e, which 
depends on the concavity of the kernel G, the solution is not fully concave on its support 
any more, but bell shaped, compare Figure [l] (c). 

To make this result more clear we present in Figure [2] the stationary solutions of the 



evolution equation (2.1) for e E (0, 1). We consider a compactly supported initial datum 
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(a) 




Figure 1. Largest eigenvalues and corresponding eigenfunctions of the 
operator Ql 4.11 (a) Largest eigenvalues e = e(L) of on L = (0, 20]; (b) 
Corresponding eigenfunctions for e(L) with L G [0,1]; (c) Corresponding 
eigenfunctions for e(L) with L G [1,8]. 



Figure 2. Stationary solutions for equation (2.1) with e G (0, 1 



0) = po with unit mass po = 1. As the results are the same as in Figure nl (b)-(c), 
we recognize again that up to a certain e the solutions are concave and then bellshaped. 
Furthermore we have mass conservation. As proven before in this paper, for e > J G = 1 
we do not have steady states because the impact of the diffusive term is higher then of 
the aggregation term. In this case we expect the solutions to behave like the self-similar 
Barenblatt-Pattle profiles. For e = we obtain an unique stationary solution (with zero 
center of mass), which is a Dirac-(^-distribution with unit mass centered at zero. 



Furthermore we computed the stationary solutions (2.1) with the kernel G{x) = 
|exp(— |x|) (see Figure km a)) with Lipschitz singularity at the point zero. We present 
the results in Figure [sj^bj for different e. Additionally we present eigenfunctions for the 
corresponding operator Ql 4.11| in Figure [sjj^c) , (d) . 
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Figure 3. (a) Kernel G{x) = ^exp(— |x|); (b) Stationary solutions for 
equation (2.1) for the kernel G. (c),(d) eigenfunctions for e(-L) with L e 



[0,1] resp. L G [1,8]. 
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